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Abstract 

Sextic oscillator in D dimensions is considered as a typical quasi-exactly solvable 
(QES) model. Usually, its QES iV— plets of bound states have to be computed using 
the nonlinear and coupled Magyari's algebraic equations. We propose and describe 
an alternative linear method which is iV— independent and works with power series 
in 1/yD. Main merit: simultaneous exact solvability (for all the N QES states) in 
the first two leading orders (the degeneracy is completely removed, the unperturbed 
spectrum is equidistant). An additional merit: All the perturbation corrections 
are given by explicit matrix formulae in integer arithmetics (there are no rounding 
errors) . 
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1 Introduction 



Sextic Hamiltonian in D dimensions 

H = -A + a|rf + 6|rf + c|rf , a = a{N) 

enters many phenomenological and methodical considerations as a "next-to-solvable" 
model [|J. In fact, among all the real polynomial interactions, only the harmonic 
and sextic models can generate an arbitrary N— plet of bound state wavefunctions 
in an elementary form. All the similar models are often called quasi-exactly solvable 
(QES, cf. §). 

Unfortunately, the close parallel between the sextic and harmonic oscillator is 
not too robust and breaks down in practical applications ||. For example, the 
Rayleigh-Schrodinger unperturbed propagator ceases to be diagonal in the sextic 
case H). Moreover, the key weakness of any QES model lies in the nonlinearity of 
its secular equation which has the polynomial form of degree N |J. Non-numerical 
determination of the sextic energies is only feasible at N < 4. Otherwise, in a sharp 
contrast to harmonic case, the values of energies E n are only available up to some 
rounding errors. 

In order to refresh the parallels we shall describe a new approach to the sextic 
QES bound state problem. It is based on some surprising results of the symbolic 
manipulation experiments. They were performed in MAPLE using the technique 
of Groebner bases. We revealed that the QES energies become equidistant and 
proportional to integers in the limit of the large spatial dimensions D — > oo. This 
feature is presented in Sections 2 and 3. 

In the second step of our analysis one discovers that the systematic evaluation of 
the Rayleigh-Schrodinger corrections proves feasible in closed form. In spite of the 
non-diagonality of propagators, a merely slightly modified form of construction can 
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be used. It gives the energy formula 

E(X) = + A + A 2 + . . . + \ K + 0(\ K+1 ), A = 1/VD. 

Its coefficients E^ are obtainable without any rounding errors (cf. Sections 4 and 5 
below). 



2 An unusual solvable limit: Large dimensions D 



All the sextic oscillator states are determined by the radial Schrodinger equation 

^(r)=E^(r). (1) 



d 2 £(£+1) 2 , 4 6 

+ v - ' +ar 2 + br 4 + cr b 



dr 

It contains the dimension D and the angular momenta k — 0, 1, . . . in £ = k + (D — 
3)/2. The elementary ansatz 

00 / 1 1 \ 

${f) = h n r 2n+e+1 exp( — /3r 2 -- 7 r 4 , c = 7 2 > 0, b = 2/3 7 > (2) 

n=0 A J 

converts this ordinary differential equation into the linear algebraic system charac- 
terized by the tridiagonal Hamiltonian matrix, 



QWh = Eh, QW 



Bo C 

A 1 B x d 



V 



(3) 



An-2 -Bjv_ 2 Cn-1 

where the dimension is to be infinite, iV — > oo, and the matrix elements are elemen- 
tary, 



A n = 7 (An + 2£ + 1) + a - (3 2 , B n = B n (E) = (3 [An + 2£ + 3), 
C n = -2(n + 1) (2n + 2£ + 3), n = 0, 1, . . . . 



(4) 



2 



The (quasi-)variational limit iV — > oo gives the numerically correct spectrum ||. 
For the sake of simplicity, let us now constrain our attention to the simplified model 
of Singh et al || characterized by the QES condition inposed upon the quadratic 
coupling a = a(N), 

a(N) 



1 

47 



2 6 2 -7(4iV + 2£+l). 



In this way one achieves the rigorous termination of the wavefunctions, 



w+i 



in +2 



0. 



(5) 



The latter assumption merely changes the lower diagonal in eqs. (|3|) and @ to the 
shorter formula A n = 47 (n — N). Exact energies become available only at the first 
few integers N < 4. Beyond N — 4, QES solutions remain numerical. Moreover, the 
intrinsic asymmetry of our Hamiltonian ([3]) causes a loss of precision which grows 
quickly with the degree N |J. 

In such a setting we have noticed, purely empirically, that the solutions are getting 
simpler when the spatial dimensions grow, D> 1. In the leading-order approxima- 
tion, the corresponding matrix Schrodinger equation becomes diagonally dominated, 



E-/3D 2D 
A(N - 1)7 E — (3D AD 



hi 



V 



67 E-/3D 2(N-1)D 
4 7 E-/3D J 



hN-2 

1 hN-i J 



0. 



(6) 



This enables us to evaluate the fully degenerate dominant eigenvalue, 



E = pD-2J2-fDz 



(7) 



where z is a constant. 
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3 The removal of degeneracy in sub-dominant ap- 
proximation 

Once we switch to the new energy variable z, we may pre-multiply eq. @ by a 
diagonal and regular matrix with elements fP where p = \J D / (27). This leads to 
the new, non-diagonal matrix Schrodinger equation. It determines the leading-order 
components of the renormalized Taylor coefficients pj = [D/(2 r y)Y' 2 hj and has the 
following transparent form, 





(N-l) 



1 





V 



(JV-l) 

1 



Po 
Pi 

PN-2 

\ Pn-i J 



Po 
Pi 

PN-2 

\ Pn-i J 



In spite of the manifest asymmetry of this equation, all its eigenvalues remain strictly 
real. We computed these eigenvalues by symbolic manipulations in integer arith- 
metics and discovered that the underlying nonlinear secular equation is solvable 
exactly and completely. The iV— plets of its energy roots proved nondegenerate, 
equidistant and extremely elementary, 



[z u z 2 , z 3 ,..., z N -!, z N ) = (-N + 1, -N + 3, -N + 5, . . . , N - 3, N - 1] 



(9) 



This result is valid at an arbitrary finite matrix size N. 

It is quite elementary to verify that also the the respective left and right eigen- 
vectors remain real. Up to their norm, all of them can be represented in terms of 
integers. Their components may be arranged in the rows and columns of certain 
square matrices, 

P(0) = 1, 
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P{2) 



v/4 



1 
2 
1 



1 1 

-2 
-1 1 



P(4) 



16 



P(3) 



1 

2 


-2 
-1 



1 


-2 

1 



1 

-2 

2 

-1 



1 1 

3 1 

3 -1 

1 -1 

-4 
6 
-4 



-3 
3 



1 -1 



etc. These matrices P = P(N — 1) are all asymmetric but idempotent, P 2 = I. 

We may summarize that in the limit D — > oo, the QES sextic model may be fac- 
torized easily. After a suitable normalization, all the components of the eigenvectors 
are integers. 



4 An adapted Rayleigh-Schrodinger perturbation 
recipe 

At the finite values of D and starting directly from the second- order precision of pre- 
ceding section, the routine perturbation theory becomes applicable since the unper- 
turbed Hamiltonian remains diagonal and all its spectrum is safely non-degenerate. 

At any » the Schrodinger equation (^) is an eigenvalue problem with the 
perturbed Hamiltonian of the two-term form, 

H{\) = P (0) + A H {1) + A 2 P (2) , A = 1/VD. 

Both the perturbations are one-diagonal matrices which depend on the value of the 
angular momentum k, 

=JL(2n + k), n = 0,l,...,A-l, 
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H {2) ) nn+i = -(n + l)(2n + 2k), n = 0, 1, . . . , N - 2. 

We may re-write our Schrodinger equation (|3|) in the textbook perturbation-series 
representation at any N, 

(ff(°> + A#« + A 2 #( 2 >) ■ (V (0) + A^ (1) + . . . + A*V (J ° + 0{\ K+l )) 
= (^<°> + . . . + A K V (X) + 0(A^ +1 )) ■ (e(°) + . . . + A*eW + (9(A X+1 )) . (10) 

Let us again concatenate the (lower-case) zero-order vectors p = = ij)^ into an 
N by N matrix P = p(°\ with all the eigenvalues arranged also in a diagonal matrix 
e^°\ In this way the zero-order equation H^tp^ = ip^e^ is satisfied identically. 
Indeed, in our compactified notation, it reads Pe^PP = Pe^ and we know that 
P 2 = L 

With the factorized # (0) = Pe^P, we shall use the same convention in all orders 
and concatenate the vectors iff- , j = 1, 2, . . . , N in the square matrix ^W. In the 
first order of perturbation analysis this replaces the O(X) part of eq. ( |10D by the 
matrix relation 

e (i) + P *(i) e (o) _ £ (0) P gr(i) = P p(D p. ( U ) 
In the second order we get 

e P) + p *(2) e (0) _ £ (0)p (2) = p p(2) p + p #(1) ^(D _ p ^(1) £ (1) (12) 

etc. The available expressions occur on the right-hand side of these equations while 
the unknown quantities stand to the left. All the higher-order formulae have the 
same structure. 

We may summarize that the diagonal part of equations (11) or (|Tj|) determines 



the energy corrections e^ 1 ' and e^ 2 \ respectively. Non-diagonal components of these 
matrix relations are to be understood as a definition of the eigenvectors. 
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5 Merits of the method: An N 



= 2 illustration 



One has to move up to the higher-order level for the elimination of the normalization 
ambiguities. This has been multiply clarified in the literature on perturbation theory 
0. Still, we should emphasize a user-friendliness of this normalization freedom 
within the framework of the present formalism. For illustration, let us consider just 
the s— wave problem in the N = 2 case. Immediately, our first-order formulae give 
the two energy corrections which are both equal to each other, 

= eg = P/yfey. (13) 

One discovers that the O(X) level of precision provides just an incomplete information 
about the norms of the first-order wave functions. This is the well know normalization 
freedom manifesting itself in the present setting. On the O(X) level of precision only 
two constraints — = ~P/V^1 an d ^12 + ^22 = P/V^l are imposed upon 
the wavef unctions. Their definition must be completed in the subsequent order. 

In any higher order computation, the use of the computerized symbolic manip- 
ulations is strongly recommended. Their implementation is trivial. The algorithm 
can be written in integer mathematics and generates, therefore, the perturbation 
series without any errors. This is our most important conclusion. One generalizes 
immediately the above leading-order results (|7|), @ and (p~3|) to the compact energy 
series for our particular sextic k = N — 2 = illustration, 

E 1>2 = ^±^ + 2p±J^\ + 0.\ 2 T J^\ s + 0.\* + O(\ 5 ). (14) 
A z A ^27 87A/27 

One can observe the (complete) leading-order degeneracy of section 2 as well as 

its immediate next-order removal (|9|) as discussed in section 3. It is also amusing to 

notice the above, hand-evaluated and quite unexpected, degeneracy of the subsequent 

0(1) correction. 

One can notice the existence of certain identically vanishing corrections here. In 
fact, their rigorous evaluation would not be possible within the standard framework 
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of perturbation theory where the summations over the intermediate states must be 
computed in finite precision. Only within the present formalism which is able to work 
in integer arithmetics, the unusual feasibility of proving the precise cancellation of 
the series of corrections can be achieved. This is one of the less expected though 
most important merits of our present methodical proposal and construction. 
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